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I.  Introduction 


A  p  r  ‘  s 


In  a  previous  paper  ([£])  one  of  us '-presented  an  iterative 
method  for  solving  the  exterior  Dirichlet  problem  for  the  Helm¬ 
holtz  equation  defined  in  the  plane  and  used  this  result  to  pro¬ 
vide  a  constructive  approach  for  solving  the  low  frequency  in¬ 
verse  scattering  problem  for  a  cylinder.  These  results  were 
based  on  the  use  of  conformal  mapping  and  the  fact  that  the 
integral  of  the  normal  derivative  of  the  total  field  over  the  bound 
ary  of  the  obstacle  vanishes,  neither  of  which  is  valid  in  the 
three  dimensional  case.  ,  in-  this  paper,  we  shall  show  how  the  k 

analysis  'of  [53^'can  ue  modified  in  order  to  extend  these  results 
to  the  case  of  the  exterior  Dirichlet  problem  for  the  Helmholtz 
equation  in  k  .  Our  results  are  based  on  chosing  an  appropriate 
fundamental  solution  such  that  the  integral  equation  associated 
with  the  exterior  Dirichlet  problem  can  be  solved  by  iteration 
for  sufficiently  small  values  of  the  wave  number. ^  Previous  results 
in  this  direction  have  been  given  by  Kleinman  ([9])  and  Ahner 
([!]).  However,  in  Kleinman’ s  approach  it  was  necessary  to  first 
construct  the  Green's  function  for  Laplace's  equation  defined 
in  the  exterior  of  the  scattering  obstacle  D,  whereas  in 
Ahner 's  approach  it  was  necessary  to  compute  the  first  eigen¬ 
function  of  the  integral  equation  associated  with  this  problem. 

Our  method  avoids  both  of  these  computations  and  instead  provides 
an  integral  equation  formulation  of  the  exterior  Dirichlet  problem 
for  the  Helmholtz  equation  with  the  kernel  of  the  integral  equa¬ 
tion  being  independent  of  D  and  expressable  in  closed  form.  Such 
a  formulation  is  particularly  suitable  for  obtaining  analytic 


approximation  for  low  frequency  scattering  problems  (c.f.  [2]). 
However,  rather  than  presenting  an  application  to  the  direct 
scattering  problem,  we  shall  show  how  our  iterative  procedure 
can  be  used  to  obtain  information  on  the  inverse  scattering  problem 
for  low  frequency  acoustic  waves .  The  importance  of  our  itera¬ 
tive  procedure  in  this  case  is  that  it  allows  us  to  rigorously 
establish  the  limiting  behaviour  of  the  total  field  as  the  wave 
number  k  tends  to  zero.  This  step  is  in  general  nontrivial 
since  in  the  classical  formulation  of  the  scattering  problem  as 
an  integral  equation  the  resolvant  operator  has  a  pole  at  k  =  0 
and  we  are  thus  faced  with  a  singular  perturbation  problem  (c.f. 

[12]  and  the  references  contained  therein.  In  this  regard  we 
note  that  in  many  cases  it  is  not  in  fact  the  limiting  behaviour 
of  the  total  field  that  is  of  interest,  but  instead  the  coef¬ 
ficients  of  certain  higher  order  terms  -  c.f.  [4],  [5]).  Having 
established  the  limiting  static  problem  associated  with  the  in¬ 
verse  scattering  problem,  we  proceed  to  follow  a  standard  approach 
(c.f.  [16])  in  order  to  obtain  information  on  the  shape  of  the 
scattering  obstacle  from  a  knowledge  of  the  low  frequency  behaviour 
of  the  far  field  pattern.  Our  contribution  here  is  to  note  that 
although  in  general  the  far  field  pattern  is  only  known  approxi¬ 
mately,  the  determination  of  the  shape  of  the  obstacle  D  is 
stabilized  by  the  fact  that  the  solution  of  the  static  potential  pro¬ 
blem  is  known  to  have  values  lying  between  zero  and  one,  i.e.  the 
potential  function  to  be  determined  lies  in  a  compact  set  of  har¬ 
monic  functions  (c.f.  [8],  [15]).  Here  it  is  necessary  to  assume 
that  we  know  "a  priori"  the  radius  of  a  ball  containing  D  in 


3 


its  interior,  with  optimal  results  holding  if  we  know  the  radius 
of  the  smallest  such  ball  containing  D. 


II .  The  Exterior  Dirichlet  Problem  for  the  Helmholtz  Equation 
We  begin  by  first  considering  the  exterior  Dirichlet  problem 
for  Laplace's  equation,  i.e.  to  find  u  «  C2(R3\D)  n  C°(13\D) 
such  that 


(x)  +  us(x)  in  S3\D 

(2.1a) 

u  s  0  in  R3\D 

(2.1b) 

=  0  for  x  c  3D 

(2.1c) 

lim  us(x)  =  0 

(2. Id) 

•  .  3 

where  D  is  a  bounded  simply  connected  domain  in  R  containing 

.  2  i 

the  origin  with  C  boundary  3D,  u  is  a  given  solution  of 

Laplace's  equation  defined  in  all  of  R3,  and  r  *  |x| .  We  note 

•m 

that  under  these  conditions  a  solution  u  to  (2.1)  exists  and 
u  e  C2(R3\I))  n  C1(R3\D)  (c.f.  C7]).  Now  let  B  be  a  ball  of  radius 
R  contained  in  D  with  center  at  the  origin  and  let 


G(x,y) 
■»»  ** 


1  R  1 

|x-y  I  "  TxT  [xA-y I 


(2.2) 


be  the  Green's  function  for  Laplace's  equation  in  the  exterior  of 
B.  Then  from  Green's  formula  we  have 

4?  j  TvTyT  “  uS  JvTyT  S(*»y)}  da(y) 

3D 

(2.3) 

-1/2  us(x),  x  «  3D 

8  ", 

— u™ ( X ) , X  €  RJ\P 
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where  v  denotes  the  unit  outward  normal  to  3D.  In  addition, 
since  u1  is  a  solution  of  Laplace’s  equation  in  all  of  R3 
and  G(x,y)  =  0  for  y  e  33,  we  have 

w  I  {«*■!!>  w  - ^  w  s'?’?'} 

3D 


+ 


G(^’y)  ds(y) 


' 

1/2  u^Cx) ,  x  £  3D 
0  ,  X  £  R3\D 


(2.4) 


where  again  v  denotes  the  unit  outward  normal.  Adding  (2.3) 
and  (2.4)  together  and  using  the  fact  that  u(x)  =  0  for  x  £  3D 
now  gives 


-u1(x) 


!  s(i’V  3v7yT  is(y>  *  5?  I 

3D  33 


u 


JvfyT 


G(x,y)  ds(y) 


3  __ 

-u(x)  ,  x  £  X  \D 
-1/ 2u(x) ,  x  £  3D 


From  the  jump  discontinuity  properties  of  the  single  layer  potential 
we  now  have  the  following  integral  equation  for  3u/3v  evaluated 
on  3D: 


3v 


+  2?  I  Tuhr  1^777  ds(y) 

3D 


2  377S7  -  FvfxT  77T  I  ul  TvTyT  G(?’y)  ds(y)  }  *  e 

3B 

N- 


(2.6) 
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Not*  that  in  the  special  case  where  u1(x)  =  1  we  have 

Iff  |  TvTyT  G(^^)  ds(y)  s  2  tIt  (2*?) 

3B 

for  x  c  3D  and  hence  in  this  case  the  right  hand  side  of  (2.6) 
becomes  -2  R  jJjJy  . 

We  shall  now  show  that  the  integral  equation  (2.S)  cam  be 
solved  by  successive  approximations.  (The  solution  of  (2.1) 
can  then  be  obtained  by  substituting  the  solution  of  (2.6)  into 
(2.5),  i.e.  we  have  a  constructive  method  for  solving  (2.1).) 

To  this  end  let  COD)  denote  the  Banach  space  of  continuous 
complex  valued  functions  defined  on  3D  and  equipped  with  the 
maximum  norm  and  define  the  compact  integral  operator 
K  :  COD)  -  COD)  by 

( Kv )  ( x )  a  -  ^  |  G(x,y)  v(y)  ds(y)  .  (2.8) 

3D 

If  we  let  f(x)  be  the  right  hand  side  of  (2.6)  then  (2.6)  can 
be  written  as 

v  -  Kv  a  f  (2.9) 

•a 

where  v  a  3u/3v.  In  order  to  show  that  (2.9)  can  be  solved  by 
iteration  w*  shall  show  that  the  spectrum  of  K  is  contained  in 
(•1,1).  We  first  note  that  since  K  is  compact  all  spectral 
values  of  K  different  from  zero  are  eigenvalues.  Let  X  be  an 
eigenvalue  of  K  corresponding  to  the  eigenfunction  $ , 
a  ,  and  define 
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u(x)  =  T?  j  3(x»v).$(y)  ds(y)  ;  x  £  5  \B 


(2.10) 


2 _ _ 

Then  u  is  harmonic  in  R  \D  and  DXS-,  and  from  the  jump 
discontinuity  properties  of  single  layer  potentials  we  have 


u+  =  u  on  3D 


5u+ 

3TT  5  '  £  $  +  * 
=  (-X  +  l)(j> 


(2.11a) 


(2.11b) 


on  3D 


where  +_  denote  the  limits  from  outside  and  inside  D 
respectively.  Hence 


3u . 


3u 


(1.X)  _  +  (1+X)  ^  =  0 


(2.12) 


and  since  u(x)  s  0  for  x  c  3B  we  have  from  (2.11a),  (2.12) 
and  Green's  theorem  that 


0  =  (1-X) 


'  -  3u+  r  - 

3u  C 

u+  JV~  ds  +  <1+x)  j  u- 

ds  -  (1+X) 

3D 

3B 

=  (X-l) 


If' 

R3\D 


[grad  u| ^dx  +  (X+l) 


grad  u |  dx  . 


DVB 


(2.13) 


Now  define 


I; 


I(u)  =  j  j  | grad  u | ‘dx 
D/B 


I( 


-‘If  ' 


grad  u|  dx 


R3\D 


(2.14) 
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Then  if  I(u)  +  Ku)  >0  we  have 


X  =  ~  Ku) 

i(u)  +  I(u) 


(2.15) 


which  implies  that  in  this  case  all  the  eigenvalues  of  K  are  real 
and  contained  in  (-1,1).  If  x(u)  +  I(u)  =  0  then  u  is  constant 
in  S^\5’  and  D\3  and  since  from  (2.10)  we  have  that  u  vanishes 
at  infinity  and  on  33  we  can  conclude  that  u(x)  =0  for 
x  s  R^\3.  Hence,  from  (2.11b)  we  have 

f 3u  3u, 


„  1 
‘  2 


-  0 


(2  .15) 


i.e.  X  is  not  an  eigenvalue,  contrary  to  assumption.  Hence, 
the  spectrum  of  K  is  contained  in  (-1,1)  and  hence  (2.9)  cam 
be  solved  by  successive  approximations. 

Before  turning  our  attention  to  the  exterior  Dirichlet  pro¬ 
blem  for  the  Helmholtz  equation  we  briefly  consider  what  effect 
a  change  of  the  radius  R  of  -the  ball  3  has  on  the  spectral 
radius  of  the  operator  K.  In  particular,  consider  the  special 
case  when  D  is  a  ball  of  radius  one  and  u  is  defined  by  (2.10). 
Then  u  has  an  expansion  in  spherical  harmonics  of  the  form 

Z  E  (anm  pn  +  b^,  p"”"1)  (cos  e)eim*;  R<p<l 
n=0  m*-n 


u(x)  = 


•  n 

Z  Z 
n=Q  m=-n 


nm  "  ~nm 
-n-1  ^m 


c  0 
nm  w 


(cos  6)e^jn^;  l<p<» 


(2.17) 


where  (p,9,$)  are  spherical  coordinates,  P®  denotes  an  as¬ 
sociated  Legendre  polynomial,  and  the  a _ ,  b_  .  c  _  are  con- 

nm  nm  nm 

stants.  Then  (2.11a),  (2.12)  and  the  fact  that  u(x)  »  0  for 
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and  this  system  has  a  nontrivial  solution  if  and  only  if  the 
determinant  of  the  coefficients  vanishes,  i.e. 


_  l-( 2n+2)R2n*1 
n  "  2n+l 


C2.19) 


From  (2.19)  we  see  that  in  this  special  case  the  eigenvalues 
decrease  as  R  increases,  i.e.  by  using  overrelaxation  we  can  make 
the  spectral  radius  smaller.  This  in  turn  implies  that  the  rate 
of  convergence  of  the  successive  approximations  to  the  solution 
of  the  integral  equation  is  improved.  By  using  the  methods  intro¬ 
duced  in  [13]  it  can  in  fact  be  shown  that  this  behaviour  is  true 
not  only  for  a  sphere  but  also  for  an  arbitrary  bounded  simply 
connected  domain  ([11],  [19]).  Hence  for  computational  purpose 
it  is  desirable  to  pick  R  as  large  as  possible  and  then  use  over- 
relaxation  to  improve  the  rate  of  convergence  (c.f.  [13]). 

We  now  consider  the  exterior  Dirichlet  problem  for  the 

2  3  —  3 

Helmholtz  equation,  i.e.  to  find  u  e  C  (R  \D)  n  C°(R  \D)  such 


that 
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u(x)  =  ui(x)  +  us(x) 

in  R3\D 

(2.20a) 

n 

A3u.  +  k“u  =  0  in 

r3\S 

(2.20b) 

u(x)  =  0  for  x  ( 

•  3D 

(2.20c) 

.  .  f 3uS  s I 

lim  r  v—  -oku 

r-wo  l3r  J 

=  0 

( 2 . 20d) 

where  the  "incoming  wave"  u1  is  a  solution  of  (2.20b)  in  all  of 
R3,  the  wave  number  k  is  positive,  and  the  radiation  condition 
(2.20d)  for  the  "scattered  wave"  us  is  assumed  to  hold  uni¬ 
formly  in  all  directions.  The  domain  D  is  assumed  to  satisfy 
the  same  conditions  as  in  problem  (2.1)  and  under  these  condi¬ 
tions  we  can  again  conclude  that  a  solution  to  (2.20)  exists  and 
u  £  C2(S3\D)  n  C1(S3\D)  ([18]).  Now  define 


GkCx,y) 


eik!x-x| 

fx-yr" 


ik |x*-y 
Ix^'y r 


(2.21) 


where  R  is  again  the  radius  of  the  ball  3  contained  in  the 
interior  of  D.  Then  with  respect  to  y  the  function  G.  is 
a  fundamental  solution  to  the  Helmholtz  equation  satisfying  the 
radiation  condition  (2.20d).  Hence,  proceeding  in  the  same 
manner  as  we  did  for  Laplace's  equation,  we  can  use  Green's  for¬ 
mula  to  derive  the  following  integral  equation  for  the  unknown 
normal  derivative  of  the  solution  to  (2.20)  evaluated  on  3D: 


3u 

TvTxT  T  2tt 


3v 


?  \  TvTxT  VS’*’  JvTyT  ds(y)  5  2  TUTSI 
3D 

1ST  2?  I  (w  Gk(*’y)  -  ul  TvTyT  Gk(*’y)}  ds(y)’  5 


(2.22) 

3D. 


+ 


Note  that  the  additional  term  on  the  right  hand  side  of  (2.22) 

(as  compared  with  the  right  hand  side  of  (2.26))  occurs  since 
we  no  longer  have  the  fact  that  G^(x,y)  =  0  for  y  s  SB.  This 
minor  inconvenience  could  have  been  avoided  by  chosing  Gk  to 
be  the  Green's  function  for  the  Helmholtz  equation  in  the  exterior 
of  B  instead  of  as  in  (2.21),  but  in  this  case  the  kernel  of 
the  integral  equation  (2.22)  could  no  longer  be  expressed  in  a 
simple  closed  form  such  as  (2.21). 

We  shall  now  show  that  the  integral  equation  (2.22)  can  be 
solved  by  iteration  for  k  sufficiently  small.  As  in  the  case 
of  Laplace's  equation  this  provides  a  constructive  method  for  sol¬ 
ving  the  boundary  value  problem  (2.20).  Let  COD)  again  denote 
the  Banach  space  of  continuous  complex  valued  functions  defined 
on  3D  equipped  with  the  maximum  norm  and  define  the  compact 
integral  operator  :  COD)  -*•  COD)  by 

(Kkv)(x)  =  -  |  ^TTxT  Gk(*’y)  v(y)  ds(y)  *  (1*23) 

3D 

Let  f(x)  denote  the  right  hand  side  of  (2.22).  Then  the  integral 
equation  (2.22)  can  be  written  in  the  form 

v  -  Kkv  =  f  (2.24) 

where  v  =  Then  since 

|Gk(x,y)  -  G(x,y)|  =  0(k)  (2.25) 


it  follows  that 
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!  ! Kk-K I  I  s  0(k)  (2.26) 

where  |j*||  denotes  the  maximum  operator  norm.  Hence,  since  the 
spectral  radius  of  K  is  lass  than  one,  we  can  follow  the  argu¬ 
ment  of  Theorem  2.2  of  CIO]  to  conclude  that  there  exists  a 
positive  constant  kg  such  that  the  spectral  radius  of 
is  less  than  one  for  0  _<  k  <_  kQ  .  Thus  for  such  values  of  k  we 
can  solve  (2.22)  by  the  method  of  successive  approximations. 


Ill .  Remarks  on  the  Inverse  Scattering  Problem 


We  now  consider  the  application  of  the  results  of  Section 
Two  to  the  inverse  scattering  problem  associated  with  (2.20).  In 
particular,  if  we  apply  Green's  formula  to  the  solution  u  of  (2.20) 
and  the  fundamental  solution  for  the  Helmholtz  equation,  we  have 
that  for  x  «  X3\D  (c.f.  C6]) 


us(x)  = 


-  f 

4ir  J 
3D 


3u 

3v(y ) 


,ik | x-y | 
|x-y| 


,ikr 

4ttt 


F(8,0;k)  +  0 


1 


ds(y) 


(3.1) 


where  (r,9,$)  are  spherical  coordinates  and 


F(0  ,<P  ;k) 


1  (  3u  .-ikri'y 

J  TvTyT  ~ 

3D 


ds (y) ; 


n 


(3.2) 


The  inverse  scattering  problem  we  wish  to  consider  is  to  determine 
the  shape  of  3D  from  a  knowledge  of  the  far  field  pattern  F 
for  all  angles  0,$  and  low  values  of  the  wave  number  k. 
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i 

We  first  want  to  reduce  the  solution  of  the  inverse  scattering 

problem  to  a  generalized  moment  problem  for  the  conductor  potential 

ikxT 

ox  the  domain  D.  To  this  end  suppose  u~(x)  =  e  where 

x  =  (x1,x2>x3).  Then  from  the  results  of  Section  Two,  we  have  that 

the  integral  equation  (2.24)  can  be  solved  by  successive  approxi- 

/ 

mations,  i.e. 


_3u_ 

3v(x) 


(3.2) 


From  (2.26)  and  [17],  p.  164,  we  now  have  from  (3.2)  that 


3u 

avTxT 


r 

n=0 


>5  E2R 


3  1 

^  T*I 


+  0  (k) 


3u0 

3 vTO-  *  0(lc) 


(3.3) 


for  x  e  3D.  From  the  results  of  Section  Two  we  can  now  identify 

3 

uQ(x)  for  x  e  R  \D  as  the  solution  of  the  Dirichlet  problem 


uQ(x)  =  1 

-  u®(x)  in  £3\D 

(3.4a) 

A2uq  =  0 

in  £3\D 

(3.4b) 

un(x)  =  0 

for  x  c  3D 

(3.4c) 

lim  Uq(x) 

|  x  |  •+« 

=  o  , 

(3.4d) 

s 

i.e.  uQ  is  the  conductor  potential  of  the  surface  3D  (c.f. 
[3]).  Now  expand  F  in  a  series  of  spherical  harmonics 

=  Z  Z  a  (k)Pm(cos  9)  eim<^ 

_  nm  n 


F(  9  ,*}>  ;k) 


(3.5) 
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we  have  from  (3.3),  (3.9)  and  (3.10)  that  for  n  >_  0 , 
-n  <  .u  <  n. 


a  (k)  = 
nm 


k 


T(n+3/2) 


T7  j  JTTyT  £  ?S<cos  e')«-“”>]jd»<y>  *  0(k 


(3.11) 


Hence,  if  we  define 


„  -  |2n+1  r(n+3/2)  in  a  (k) 

u__  -  lim  run 

nm  k^Q  - - - 

/¥  kn 


(3.12) 


then  the  inverse  scattering  problem  has  been  reduced  to  determining 
D  from  the  "moment"  problem 

“nm  *  *  j  KI5T  E"  P"  <CCS  e,>  dSty)  •  l3-13) 

3D 

Note  that  although  the  y^  are  explicitly  computable  from  the  far 

field  pattern  F,  small  errors  in  measuring  the  coefficients 

a^  (k)  will  result  in  large  errors  in  the  numbers  y  „  if  n 
nm  nm 

is  large.  Hence,  in  practice  we  must  assume  that  only  a  finite 

number  of  the  y„„  are  known.  Observe  also  that  the  "moment" 
nm 

problem  (3.13)  is  nonlinear  since  both  uQ  and  the  region  of 
integration  depend  on  D. 

We  shall  now  show  that  (3.13)  allows  us  to  compute  uQ 
outside  of  a  ball  §  containing  D  in  its  interior.  Since  only 
a  finite  number  of  the  y^  are  known  we  can  only  compute  u^ 
approximately,  and  in  order  to  obtain  error  estimates  it  is  neces¬ 
sary  to  know  the  radius  of  B  "a  priori" .  Having  obtained  an 


15 


approximation  to  uQ ,  we  note  that  the  level  curve  uQ(x)  *  3 
is  3D  and  hence  the  level  curves  u^tx)  s  y,  0  <  y  <  1, 
approximate  3D  as  y  tends  to  zero .  Thus ,  if  we  determine 
the  curves  uQ(x)=  Y,  0  <  Y  <  1,  for  x  «  1  \§  and  extrapolate 
as  y  tends  to  zero  we  can  determine  an  approximation  to  3D. 
Obviously  the  success  of  this  procedure  i3  optimized  if  3  is 
in  fact  the  smallest  ball  containing  D  in  its  interior.  To 
determine  an  approximation  to  u^(x)  for  x  *  X3\3  we  use  the 
fact  that  uQ(x)  =  0  for  x  «  3D  and  hence  from  (3.13)  and 
Green's  formula  we  have 


’  {“0  ^  PS  <c°s  9'>  e'104’] 

A 

33 

TT  Cp”  ?n  (cos  8')  e’im<^] j-  ds. 


(3.14) 


Since  on  33, 


uQ  has  an  expansion  of  the  form 

s  l  +  E  Z  b__  o~n~l  P®  (cos  9)  e1®* 
nm  n 


(3. IS) 


we  have  from  (3.14),  (3.15)  and  the  orthogonality  of  the  associated 
Legendre  polynomials  that 


nm 


*  4t  b 


(n*m) ! 
nm  ( n-m) ! 


(3.16) 


Suppose  now  that  we  know  the 

for  -n  _<  m  _<  n,  0  <_  n  <_  N. 

N  N  n 

u”(x)  si*  I  Z 

'J  ~  n*0  m*-n 


(and  hence  from  (3.16)  the  b _ ) 

nm  nm 

We  then  want  to  know  how  accurately 

b  p~n” 1  P®  (cos  8)  e1®* 
nm  n 


•■■>•*** 


(3.17) 


approximates  u.(x)  for  x  *  JT  B.  To  this  end  we  note  that  from 


(3.15)  we  have  that  for  n  >  1  and  o  >  a  =  radius  of  3, 


nit 

u 


4*  (n+m)!  „-n-l 

H7 1  Zn^mTT  Dnm  0 


n(x)  P®  (cos  6)  e“^m^  sin  g  ded*> 
u  n 


(3.18) 


and  since  from  the  maximum  principle  0  <  uQ(x)  <_  1  for 

3 

x  «  Jt  \D,  we  have  from  (3.18)  that 


jtSett  a'n'1|l  i  2,1  C  9>l  si"  9 

1/2  _  i/: 

S  2n  I  sin  0  dfl  J  |?™(c os  d)j2  sin  9  d9 


*  23/^h 


4ir  (n+m)! 


(3.19) 


lb  j  <  a 

'  nm1  —  / 


1*1)  (n-m) ; 
[n+mj!  _ 


(3.20) 


Hence,  for  P  >  a  we  have 


1  f  iu  x 1 2 .  _  *  ?  Ih  I2  (n+m)Ip’2n~2 

~ 7  J  ds  f  ,bnm(  TrwnTTnn+IT 

4ffp  ixjsp  n=N+l  m=-n 


(3.21) 


<  2ir  E  (2n+l) 
n=N+l 


and  in  view  of  (3.16)  this  tells  us  how  many  Fourier  coefficients 
of  the  far  field  pattern  are  needed  in  order  to  approximate 
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Ug(x)  for  |x|  _>  >  a  to  obtain  a  given  degree  of  accuracy  in 

the  L"  sense.  We  note  that  considerably  sharper  results  are 
obtainable  in  the  case  of  the  inverse  scattering  problem  for  a 
cylinder  due  to  the  fact  that  conformal  mapping  techniques  are  now 
available  (c.f.  [4],  [5]). 
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